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jy^ . We extend to Lipschitz continuous functional either of the true paths or of the 

Euler scheme with decreasing step of a wide class of Brownian ergodic diffusions, the 
Central Limit Theorems formally established for their marginal empirical measure of 
these processes (which is classical for the diffusions and more recent as concerns their 
discretization schemes). We illustrate our results by simulations in connection with 
Oh ' barrier option pricing. 
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1 Introduction 



In a recent paper ( [22] ) , we investigated weighted empirical measures based on some Euler 
schemes with decreasing step in order to approximate recursively the distribution ¥ u of a 
stationary Feller Markov process X := (Xt)t>o with invariant distribution v. To be precise, 
let (X t ) be such an Euler scheme, let (Tk)k>i denote its sequence of discretization times and 
let (^/OfcM be a sequence of weights. On the one hand we showed under some Lyapunov- 
type mean-reverting assumptions on the coefficients of the SDE and some conditions on 
the steps and on the weights that 

= 1 f>F(Xr fc+ .) ^±^P,(F) = / E[F(X*Mdx) a.s., 
r)l + • • • + Vn ^ J 

7 (L1) 

for a broad class of functionals F including bounded continuous functionals for the Sko- 
rokhod topology. On the other hand, in the marginal case, i.e. when F(a) = f(a(0)), then 
the procedure converges to v{f). When the Poisson equation related to the infinitesimal 
generator has a solution, this convergence is ruled by a Central Limit Theorem (CLT): 
this has been extensively investigated in the literature (for continuous Markov processes, 
see [3], for the Euler scheme with decreasing step of Brownian diffusions, see |18t I21j). As 
concerns Levy driven SDEs, see [24] . 

Our aim in this paper is to extend some of these rate results to functionals of the path 
process and its associated Euler scheme with decreasing step, i.e. to study the rate of 
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convergence to F U (F) of (\ f*F(X s+ .)ds) t >i and (D^(oj,F)) n >i respectively. Here, we 
choose to assume that (Xt) is an M. d - valued process solution to 

dX t = b(X t )dt + a(X t )dW t , (1.2) 

where (Wt)t>o is a g-dimensional Brownian motion and b and a are Lipschitz continuous 
functions with values in M d and respectively, where q denotes the set of d x q- 
matrices. Under these assumptions, strong existence and uniqueness hold and (Xt) is 
a Markov process whose semi-group is denoted by (Pt). We also assume that (Xt) has 
a unique invariant distribution v and we denote by P„, the distribution of (Xt) when 
stationary. 

Let us now focus on the discretization of (Xt). We are going to introduce some 
continuous-time Euler schemes with decreasing step: denoting by (r n ) n >i the increasing 
sequence of discretization times starting from Tq = 0, we assume that the step sequence 
defined by j n :=T n — r n _i, n > 1, is nonincreasing and satisfies 

n 

lim 7„ = 0, r n = V^7fc U ^ °°> +oo, (1.3) 
k=l 

First, we introduce the discrete time constant Euler scheme (Xr n ) n >o recursively defined 
at the discretization times T n by Xq = xq and 

X Fn+1 = X Vn + ln+1 b(X Tn ) + a(X Tn+1 )(W Tn+1 - W Vn ). (1.4) 

There are several ways to extend this definition into a continuous time process. The 
simplest one is the stepwise constant Euler scheme (X t )t>o defined by 

VneN, vt g [r n ,r n+1 ), x t = x Tn . 

The stepwise constant Euler scheme is a right continuous left limited process (referred 
as cadlag throughout the paper, following the French acronym). This scheme is easy to 
simulate provided one is able to compute the functions b and a at a reasonable cost. One 
could also introduce the linearly interpolated process built on (^r„)n>o but except the 
fact that it is a continuous process, it has no specific virtue in term of simulability or 
convergence rate. 

The second possibility to extend the discrete time Euler scheme is what we will call 
the genuine Euler scheme, denoted from now on by (£t)t>o- It is defined by interpolating 
the two part of the discrete time scheme in its own scale (time, Brownian motion). It is 
defined by 

VnGN, We [r n ,r n+1 ), ^ t = X Tn + (t-T n )b(X Tn )+a(X Vn )(W t -W Tn+1 ) (1.5) 

Such an approximation looks more accurate than the former one, especially in a func- 
tional setting, as it has been emphasized - in a constant step framework - in the literature 
on several problems related to the Monte Carlo estimation of (a.s. continuous) function- 
al of a diffusion (with a finite horizon) (see e.g. [7], Chapter 5). This follows from the 
classical fact that the L p -convergence rate of this scheme for the sup norm is y^y instead 
of — y / 7log7 for its stepwise constant counterpart (where 7 stand for the step). On the 
other hand, the simulation of a functional of (Ct)te[r,r+Tl is deeply connected with the 
simulation of the Brownian bridge so that it is only possible for specific functionals (like 
running maxima, etc). 
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A convenient and synthetic form for the genuine Euler scheme is to write it as an Ito 
process satisfying the following pseudo-diffusion equation 



tt = x Q + f b(Qds+ f a(QdW s (1.6) 
Jo Jo 

where 

t = T N{t) with N(t) = min{n > 0, r n+i > t}. (1.7) 

Taking advantage of this notation for the stepwise constant Euler scheme, one can also 
note that 

VieR+, X t = X t _. 

When necessary, we will adopt the more precise notation I 1 ''' 1 "' for a stepwise constant 
continuous-time Euler scheme to specify starting at x E R rf at time with a nonincreasing 
step sequence (h n ) n >i satisfying (jl.3D . 

Since we will deal with possibly cadlag approximations of continuous processes we will 
introduce the spaces ]D uc (_r, R ) of R^-valued cadlag functions on / — R+ or [0, T 1 ], T > 0, 
endowed with the topology of the uniform convergence on compact sets, rather than the 
classical Skorokhod topology (see |6J). In fact, one must keep in mind that if a : / — > R rf 

is a continuous function and (a n ) is a sequence of cadlag functions, a n — )■ a iff a ^ a 
(with obvious notations). Furthermore, usual Skorokhod distance d$k (so-called J\ and 
J 2 topologies) on D([0, T], R d ) all satisfy 

dsk(<*,P) < \\<x-P\\ T ■■= sup \a(t)-/3(t)\ 

te[o,T] 

so that any functional F : D([0, T],IR d ) — > R which is Lipschitz with respect to such a 
distance dsk will be Lipschitz continuous with respect to || . || T (hence measurable with 
respect to the Borel u-field induced by the Skorokhod topology). 

At this stage, we need to introduce further notations related to the long run behaviour of 
processes (or simply functions). Let 6 a (d/3) denotes the Dirac mass at aG lD(R + ,R rf ) and 
c/ u ) := (a u+ t)t>o denotes the it-shift of a. 

We will see below that our aim is to elucidate the asymptotic ¥(duj)-a.s. weak behaviour 

If 1 

of the empirical measures - / 6 Y (s)^(d(3)ds as t goes to infinity, where Y will be the 

t Jo 

diffusion X itself or one of its (simulatable) Euler time discretizations. This suggests to 
introduce a time dicretization at times T n of the above time integral like we did to define 
the Euler scheme. This leads us to introduce, for any aG D(R + ,R d ), the following abstract 
"Euler" empirical means 

1 n 1 f r " 

n k=l L n JO 

Then, for a functional F defined on B(R + ,R d ) and a € B(R+,R d ), 

r i n 1 fF n 

v {n) (a,F) = / F(J3)d n \a,dp) = — V lk F{a^-^) = — / F(a&)da. 

Jn(R + ,R d ) n k=i n 

In the following, we will use this sequence of empirical measures for both stepwise con- 
stant and genuine Euler schemes. Compared to [22J, this means that we assume that the 
sequence of weights (%) satisfies r\ n = j n for every n > 1. 
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Additional notations. D> (x,y) = ^2i x iVi wn l denote the canonical inner product and 
\x\ = \J (x, x) will denote Euclidean norm of a vector xG R d . 

t> Let A = [a™] G g be an R-valued matrix with d rows and q columns. A* will denote 

the transpose of A, Tr(A) = J2i a a its trace and \\A\\ := y/Tr(AA*) = {J^ijtfj)^- If 
d = q, one writes Ax® 2 for x*Ax. 

2 Main results 

2.1 Assumptions and background 

We denote by (J r t)t>o the usual augmentation of a(W s ,0 < s < t) by P-negligible sets. Since 
b and a are Lipschitz continuous functions, Equation (|1.2|) admits a unique (J r t)-adapted 
solution [Xf)t>o starting from x € R . More generally, for every u > and every finite 
^-measurable random variable H, we can consider (Xj '^) t > , unique strong solution to 
the 5L>S: 

dy 4 = b(Y t )dt + (i(y t )^ (lt) , y = 3, (2.8) 

where = W u+ t — W u , t>0, is the li-shifted Brownian motion (independent of Fu). 

Note that Xf = xf ] ' x and that X t (u) ' S can be also defined through the flow of (jl.2p by 
setting 

X (u),B ={x (u U) ^ 

Throughout this paper, we consider a measurable functional F : D uc ([0, T], M d ) —¥ M.. We 
will denote by F T the stopped functional defined on 1D) UC (M + , M. d ) by 

Va G D uc (IR + ,M a! ), F T (a) = J F(a T ) with a T (t) = a(t A T), t>0. (2.9) 

Let us introduce the assumptions on F. 

(Cp): F : D ttc ([0,T],R a! ) — > M. is a bounded and Lipschitz continuous functional. 
We set 

f F (x) = E[F T (X X )} = E[F(X? ,0 < t < T)]. 

It is classical background (see e.g. [15]) that, under the Lipschitz assumption on b and a, 
E[sup tg [ 0j T] \Xf — Xf\] < Cb l<T ,T\x — y\ so that f F is in turn clearly Lipschitz continuous. 
Additional regularity properties (like differentiability) can be transfered from fp provided 
F, b and a are themselves differentiable enough (see e.g. [IS]). Furthermore, it follows 
from its very definition and the Markov property that 

u(f F ) = F V (F T ) = J E[F T (X x )}v(dx). 

(Cp): There exists a bounded C 2 -function g F : M. d — > R with bounded Lipschitz continuous 
derivatives such that 

VxGR d , f F (x)-u(f F )=Ag F 

where A denotes the infinitesimal generator of the diffusion (11.21) defined for every C 2 - 
function / on R d by 

Af(x) = (Vf,b)(x) + ~Tr(a*D 2 fa(x)). 
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Remark 2.1. In fact, we need in the sequel that fp satisfies a CUT for the marginal 
occupation measures which follows (see |X8(, I24j) from Assumption (C F ) combined with a 
Lyapunov stability assumption (such as (S a , p ) introduced below). Namely, we have for a 
class of regular functions / satisfying / = Ag + C 



f{Xf)ds - — ^ N (0, a}) (2.10) 



, \ 'Ju n—t+oo „ 

and as soon as > —== > 0, 

(f-tl.nxr^) -,(/)) -£->JV(0,.J), (2.11) 



where 



/ |a*V5(x)| 2 ^(dx) = -2 / g(x)^(x)^(dx) 



and C denotes the weak convergence of (real valued) random variables. For details on 
results in these directions, see [3] for the continuous case and [THJ [2H E3] f° r the decreasing 
step Euler scheme. 

Checking when Assumption (C F ) is fulfilled is equivalent to solve the Poisson equation 
Au = f on When / has compact support, well-known results about the same equation 
in a bounded domain lead to Assumption (C F ) when the diffusion is uniformly elliptic (see 
e.g. [16J, Theorems III. 1.1 and III. 1.2). Such an assumption on f F is clearly unrealistic. 
In the general case, in [26], [27] and [28], the problem is solved under some ellipticity 
conditions in some Sobolev spaces and controls of the growth are given for u and its first 
derivatives. Finally, when the diffusion is an Ornstein-Uhlenbeck process, one can refer to 
|18j where the problem is solved in C 2/ 



Let us now introduce the Lyapunov-type stability assumptions on SDE fjl .2f) . Let £Q(R d ) 
denote the set of Essentially Quadratic functions, that is C 2 -functions V : K d — > (0, oo) 
such that 

lim V(x) = +oo, |W| < CW and D 2 V is bounded. 

\x\— >+oo 

Note that since V is continuous, V attains its positive minimum v > so that, for any 
A, r > 0, there exists a real constant C A r such that A + V r < C A r V r . 
Let us come to the mean-reverting assumption itself. First, for any symmetric d x d matrix 
S, set := max(0, Ai, . . . , A^) where Ai, . . . , A^ denote the eigenvalues of S. Let a S (0, 1] 
and pE [1, +oo). We introduce the following mean-reverting assumption with intensity a: 

(S a) p) : There exists a function V G £Q(M. d ) such that: 
(i) 3C a > such that \b\ 2 + Tr(aa*) < C a V a . 

(ii) There exist E R and p > such that (W, b) + X p Tr(aa*) < - pV a , 



1 

2;ar/w(*)+(^) 

tion for the diffusion (Xt)t>o- 



where X p := - sup AT^ ^ ( r _ ^vv®vv • The function V is then called a Lyapunov func- 



In Theorem 3 of [19] . it is shown that this assumption leads to an a.s. marginal weak 
convergence result to the set of invariant distributions of the diffusion. When p > 2 and 
the invariant distribution is unique, this result reads as follows. 
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Proposition 2.1. Let a G (0, 1] and p > 2 such that (S aiP ) holds. Then, 

1 n 

sup — ^^^^fe-i) <+°° a.s. ( 2 - 12 ) 

n>l n , , 
— fe=l 

Lei denote the unique invariant distribution of (jl,2p . Then, a.s., 

^X>/(*r fc _j^±^(/) 
in fc=i 

/or every continuous function f satisfying f(x) = o(V 2 +a ~ 1 (x)) as |x| — )• +oo. 

Remark 2.2. In the case V(x) = 1 + \x\ 2 , one checks for instance that for a given 
a G (0, 1], Assumption (S a , p ) is fulfilled for every p > 1 if Tr(aa*)(x) = o(l + |a;| 2a ) as 
|x| — > +oo and 

b(x) = -p(x)^-+T(x) where dlxl 2 ^ 1 < p(x) < C 2 \x\ 2a -\ 
\x\ 

and T satisfies for every x G M. d (T(x),x) = and |7~(x)| < C(l + |z| a ). 

As concerns the uniqueness of the invariant distribution v, we need an additional assump- 
tion related to the transition P T . Namely, we assume that: 

(S^,): v is an invariant distribution for (Pt)t>o an d the unique one for P T . 
Then, u is in particular the unique invariant distribution for {Pt)t>o- In fact, checking 
uniqueness of the invariant distribution for P T at a given time T > is a standard way 
to establish uniqueness for the whole semi-group (Pt)t>o- To this end, one may use the 
following two typical criterions: 

• Irreducibility based on ellipticity: for every x G R d , P T (x, dy) has a density (p T (x, y)) y& ^d 
w.r.t. the Lebesgue measure and Xd(dy) — a.s., p T (x,y) > for every x G M. d . 

• Asymptotic confluence: for every bounded Lipschitz continuous function /, for every 
compact subset K of M. d , 

sup \P kT f(xi) - P kT f(x2)\ (see e.g. @|2I]). 

2.2 Main results 

We are now in position to state our main results. 

Theorem 2.1. Let T > 0. Assume b and a are Lipschitz continuous functions satisfying 
(S a ,p) with an essentially quadratic Lyapunov function V : R d — > (0, +oo) and parameters 
a G (0, 1] and p > 2. Assume furthermore that V satisfies the growth assumption: 

yp+a-l ( r ) 

liminf — ^ > 0. (2.13) 

|s|->-+oo \x\ 

Assume that the uniqueness assumption (S^,) holds. Finally, assume that the step sequence 
(ln)n>i satisfies (Q73p and 

3/2 

£^=<+oo. (2.14) 

k>i * k 

Let F : B uc ([0,T],R d ) -+ R be a functional satisfying (C£) and (C|). 
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(a) Genuine Euler scheme: Then 

yfa(0<- n \Z{u>),F T )--p v {F T j) —^jV(0,al), (2.15) 



where 



al = -[ IE 



1 



F 



T 



E { A 2T I F 2t) ~ E (K I F T ) - jT <J*Vg F {XZ)dW u ) 2 u(dx)^j (2.16) 



and A% := f Q (F T (X* + ) - f F (X^)) du, t > 0, (is F t+ r- adapted). 

(b) Stepwise constant Euler scheme; furthermore, if there exists 5 > such that 

7 fe 2 



^^=<+oo, (2.17) 



r n (u( n \x(u),F T )-F u (F T )) —L-+tf(o,oZ). (2.18) 



k>i v K 

then, 

c 

n— >+oo 

Remark 2.3. By a series of computations, we can obtain other expressions for a 2, . In 
particular, we check in the Appendix A that o~ 2 F reads 

u\ = 2 [ T (l - ^)C F (v)dv - 2E u (f t (X) [ T a*Vg F (X u )dW u ) + f \a*V g F (x)\ 2 u(dx), 

J0 1 V J0 J JR d 

(2-19) 

where E y denotes the expectation under the stationary regime and C F is the covariance 
function defined by 

C F (u) = E U (F T (X U+ .) - f F (X u ))(F T (X) - f F (X ))). (2.20) 

This expression is not clearly positive but has the advantage to separate the "marginal 
part" that is represented by the last term from the "functional part" which corresponds 
to the first two ones. 

For instance, when F(a) = <p(a(0)), 4> being bounded and such that (ft — v{<ft) = Ah where 
h is a bounded C 2 -function with bounded derivatives, then f F = <ft and one observes that 

the first two terms of (|2.19p are equal to so that crp= \a*X7g F (x)\ 2 u(dx). This means 

that we retrieve the marginal CLT given by (|2.11|) (under a slightly more condition on the 
step sequence which is adapted to the more general functionals we are dealing with, thus, 
more constraining than that of the original paper; see below for more detailed comments 
on the steps conditions). 

If we now consider Ft defined F T (a) = (ft(a(T)), (ft satisfying the same assumptions as 
before, one can straightforwardly deduce from a simple change of variable that the limiting 
variance is still L d |<7*V/i(x)| 2 z/(<ix). In the appendix (Part B), we show that retrieving 
this limiting variance using (|2.16j) is possible but requires some non trivial computations. 
In particular, this calculus emphasizes the intricate nature of the structure of the functional 
variance. 

Given the form of u^ n \ it seems natural to introduce the (non- simulatable) sequence 



1 f n 

- / F T (^)du 

n JO 



which in fact appears naturally as a tool in the proof of the above theorem. 
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Theorem 2.2. Assume the assumptions of Theorem \2.lV a). Then, 



^(t / F T {^)ds-V u (F T )] —±-+M{V,o-l). (2.21) 

\t JO J n^+co 

Finally, we also state the Central Limit Theorem for the stochastic process (Xt)t>o 
itself. This result can be viewed as a (partial) extension to functionals of Bhattacharya's 
CLT established in |4J for a class of ergodic Markov processes. 

Theorem 2.3. Let T > 0. Assume b and a are Lipschitz continuous functions satisfying 
(Sap) with an essentially quadratic Lyapunov function V and parameters a G (0, 1] and 
p > 2. Assume (S^) holds. Let F : B uc ([0, T], K d ) — > R be a functional satisfying (Cp) 
and (Cp). Then, for every x G M. d , 

v'tf- f F(X^,0<u<T)ds-K(F T )) (0,o- 2 F ). ( 2 - 22 ) 

This means that our approach (averaging decreasing step schemes) induces no loss of 
weak rate of convergence with respect to that of the empirical mean of the process itself 
towards its steady regime. If we look at the problem from an algorithmic point of view, 
the situation becomes quite different. First, we will no longer discuss the recursive aspects 
as well as the possible storing problems induced by the use of decreasing steps: it has 
already been done in [22] and we showed that they can easily be encompassed in practice, 
especially for additive functionals or functions of running extrema (see e.g. simulations in 
Section [7]) . 

Our aim here is to discuss the rate of convergence in terms of complexity. It is clear 
from its design that the complexity of the algorithm grows linearly with the number of 
iterations. Thus, if 7 n oc n~ p , < p < 1, then T n ~ so that the effective rate of 

convergence as a function of the complexity is essentially proportional to n 2 . However, 
the choice of p is constrained by conditions (|2.14p or (|2.17p that are required for the control 
of the discretization error. These conditions imply that p must be taken greater than 1/2 
and lead to an "optimal" rate proportional to n^~ £ for every e > 0. This means that we 
are not able to recover the optimal rate of the marginal case that is proportional to n -1 / 3 
and obtained for p = 1/3 (see [19] for details). Indeed, in this functional framework, the 
weak discretization error is generally smaller and thus, is negligible compared to the long 
time error under a more constraining step condition f)2. 14|) instead of Yl l\l V^n < +°o 
in the marginal case). 

The paper is organized as follows. In Sections [3l U] and [5j we will focus on the proof 
of Theorem 12. 11 (a) and Theorem 12.31 about the rate of convergence of the two considered 
occupation measures of the genuine Euler scheme. Then, in Section [61 we will summarize 
the results of the previous sections and will give the main arguments of the proof of 
Theorems I2.1f q) and 12.31 Finally, Section [7] is devoted to numerical tests in a financial 
framework: the pricing of a barrier option when the underlying asset price dynamics is a 
stationary stochastic volatility model. 

3 Preliminaries 

As for the marginal rate of convergence (see [18J), the first idea is to find a good decom- 
position of the error (see Lemma f 3 . 1 [) . In particular, we have to exhibit a main martingale 
component. Here, since F depends on the trajectory of the process between and T, the 
idea is that the "good" filtration for the main martingale component is (^fcT)fc>o- That 
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is why, in the main part of the proof of these theorems, we will introduce and study the 
sequence of random probabilities (V^ n ' T \u}, da)) n >± defined by: 

i pnT i n pkT 

V^ T \u,d/3) = — 5 (x) (d(5)du = — V / 6 {x) (dp)du 
where u is a deterministic real number lying in [u, u] . 

To alleviate the notations, we will denote from now on, Qk = J~kT and = E[. | Qj~], 

k>0. 

At this stage, the reader can observe on the one hand that for a bounded functional 
F, V^ n,T \uj , F T ) is Gn+i = i r ( rt+ i)T-adapted for every n > and on the other hand that 
r p\ n ' T ) (oj,F t ) is very close to the random measures i>( n \£(uj),d/3) of Theorem 12. li b) by 
taking u = u V [u] and exactly equal to its continuous time counterpart in Theorem 12,21 if 

one sets u = u. (This fact will be made more precise in Section [6|). 

Hence, the main step of the proof of the above theorems will be to study the rate 
of convergence of the sequence OP( n ' T )(w,F T )) n > to F U (F T ) for which the main result is 
given in Section [6] (see Proposition I6.3|) . In this way, we state in this section a series of 
preliminary lemmas. In Lemma |3.1( we decompose the error between this new sequence 
(■p( n > T ) (w , F T )) n >i and the target ¥ V (F T ). In Lemma [3.21 we recall a series of results 
on the stability of diffusion processes and their genuine Euler scheme in finite horizon. 
Finally, in Lemma 13.31 we recall and extend results of |19] about the long-time behavior 
of the marginal Euler scheme. 

For every k S N, we define the (/^-measurable random variable 4> F (k) by 



0,(1) = o, </> F (k) 



/ F T (C iu) )du iffc>2. (3.23) 
Jlh-i 



where Ik = [(k — l)T,kT). Please note that <j) F (k) is J-^T-measurable. 
Lemma 3.1. For every F satisfying (Cp) and (Cp), we have 

pM(u,F T ) - P„(F T ) = ^ + + Qw+1 ' 3 

where 

4 

(M n ) n >i is a {Q n ) -martingale decomposed as follows : M n = M ra j with 

i=l 

AM M = 4> F (k) - E fe _i [</>,(&)], AM k , 2 = E k [</> F (k + 1)] - E fc _i[^(fc + 1)], 

AM M = / E k ^[F T (x'~' ~)] - f F (Z u )du, AM M = - / {Vg F {S«)MSu)dW u ) 

J Ik Jik 

and (G nj i), (@ n ,2) and (O n ,3) are (Q n )-adapted sequences defined for every n > 1, by: 



Efc-i 
k=i Jlk 



F T (?~ } )-F T (X K ~ " 



k=l JL > 

e„, 2 = V ( f Ag F {i u )du - AM M ) 
k=i V ^ ~ ' 

©n,3 = {4>n{F) — E n {4> n {F))) . 
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du, 



Proof. With our newly defined notations, we have, for every n > 1, 



nT 
k=i 



Now, for every k > 1, going twice backward through martingale increments, one checks 
that 

^(k + 1) = AM^j + AMf + E fe _!(0 F (A; + 1)). 

Then, noting that Efc_i(</> F (/c + 1)) = / Kk-i(F T (^~' > ))du, we introduce the approxima- 

Jlk 

tion term A0 n> i between the genuine Euler scheme £ and the true diffusion X so that 

F (fc + l) = AM fc 1 +1 + AAf fc 2 + Ae M + / K^F^X^^du. 
At this stage the Markov property applied to the original diffusion process yields 

E k (F T (X^ x) )) = E k (EuF T (X (jiU}1 )) = E k f F (^) = f F (Z») 
since u<u< kT. As a consequence, AM| is a true £/fc-martingale increment and 

Mk + 1) = AM^i + AMf + A9 M + AM fe 3 + / f F (^ u )du 
On the other hand f F = Ag F + F V (F), so that 

/ f F (iu)du-W v {F) = f Ag F {i u )du = A6 n>2 + AM fe 4 . 
Finally, summing up all these terms yields 

since e n+ i, 3 = Mi +1 - M\. □ 

Remark 3.4. The term nj i sums up the error resulting from the approximation of 
^L)>5« ^ .^. g g u j er sc h eme ( w fth decreasing step) The term Q n ^ is a residual 

approximation term as well: indeed, if we replace mutatis mutandis (, u by X u , Ito's formula 
implies that 

9 F (X {k+1)T )-g F (X kT )= / Ag F {X u )du+ / (Vg F (X u ), a(X u )dW u ), 
so that the resulting term would be, instead of n ,2 ; 9 f( x ™t) _ Q(^j n y 
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Lemma 3.2. Let p > and T > 0. Assume that b and a are Lipschitz continuous 
functions and that there exists £Q(R ) such that \b\ 2 + ||cr|| 2 < C^^cf) for a positive real 
constant Cb^- Then, 

(i) There exists a real constant C Pt T,b,a > 0, such that for every u > and every finite 
T u -measurable random vector H 

E[ sup <fP(xf^) | F u ] < C PyTAa( jf(E) and E[ sup 4> p {^+t) \ T v ] < C P)TM ^(^)- 
te[o,T] te[o,T] 

(ii) There exists a real constant C V) t > such that, for every u > 0, 

E[ sup \Cu+t ~ | F u ] < C P>T (1 + \Cu\ P hl u)+ v 

te[o,T] v ; 

(Hi) There exists a real constant C p > suc/i that, for every n > 0, 

n sup ie„ - e r „n^rj < C p 0i(e r j7! +1 - 

«e[r„,r n+ i) 

(iu) Let p > 2. Then, there exists C Pj t,s > suc/i i/iaf, /or every u > 0, 

E P-i 

E[ sup |f u+t - £ n+t | p | jy < C p ,t4> 2 (€u)lN( u )+r 
te[o,r] v y 

Proof. The proofs follow the lines of their classical counterpart for the constant step Euler 
scheme of a diffusion (see e.g. [TJ, Theorem B.1.4 p. 276 and the remark that follows). In 
particular, as concerns (ii), the only thing to be checked is that (£ u +t)t>o is the Euler 
scheme with decreasing step 7W of X^'^ u where the step sequence 7^ is defined by 

y[ u ' = rjv(«)+i - u, 7^ = 7AT(«)+fe, k > 2. (3.24) 

□ 

Lemma 3.3. Let p > 2 and a 6 (0, 1] such that (S ajP ) holds and assume that b and a are 
Lipschitz continuous functions. 

(i) Let g : R + — > R + be a nonincreasing function such that J* °° g(u)du < +00. Let (5 k ) be 
a nonincreasing sequence of positive numbers such that J2k>i^k < +00. Then, a.s., 

r+00 

/ ElVP+^H^M^du < +00 and V S k E[VP +a ~ l (^ k _ 1)T )} < +00. (3.25) 
Jo fc>i 



(ii) W^e /lave: 

sup - / ^2+ a - 1 (g ( f s < +00 a.s. 
t>ri * Jo 



(3.26) 



sup ■ 

n>l n 



t>I 
and 

1 n 

1 -^yf+- 1 (^ (fc _ 1)r )< + oo a.s. (3.27) 
fc=i 

In particular, the families of empirical measures (j f£ 5^ s ds) t>1 and (^ J2k=i %fc-ivr) n >l 
are a.s. tight. 

(iii) Assume (Sj.). Then, a.s., for every continuous function f such that f(x) = o(Vz +a ~ (x 
as \x\ — > +<x, 

-t rt 1 n 

/ K/) and - V /(£ {ifc -i)T) K/)- 
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Proof, (i) First, note that 

roc 

J ° „>1 

where 9 n = 7" 1 f^ n ^ g(u)du. Consequently, the first statement is simply a rewriting 
with continuous time notations of Lemma 4 of [19j. As concerns the second one, using 
Lemma l3.2( i) with <fi = V and the exponent f+ a_1 yields for every k > 1 and every u G 

Et^- 1 ^)] < Cp,a,TE[^ +a - 1 (^)]- 

As a consequence, considering the integrable, nonincreasing, nonnegative function g = 
Efc>i 1 h-i 5 k leads to 

YmV^-H^k-DT)} < C p , a ,rY I W P+a ~\Q\9{n)du < +00 



owing to the previous statement. 

(ii) Set r=|+a— 1>0 since p > 2 and a > 0. First, for every n > 1 and every 

t £ [r n , r n+ i), 



n+l 



n+l 



1 />t r 1 9 

since 7„ is nonincreasing. Now, owing to Proposition ^. 11 

1 ™ 

^ 1 J- n 



n> 



fc=l 



and (|3.26|) follows. 

Let us deal now with ()3.27|) . Given (13.26f) . it is clear that (|3.27p is equivalent to showing 
that for an increasing sequence (t^) such that to = 0, sup fc>1 (ifc — ifc_i) < +00 and 
tk ->■ +00, 



sup - V ( (t fc - ^-i)U r (£ fcT ) - / U r (£ u )du ) < +00 a.s. 



(3.28) 



Setting tfc = T N ( kT } + i for every > 1, this suggests to introduce the martingale defined 
by Nq = and for every n > 1, 



1 



k=l 



tk-l 



where E 4fc [.] := E[ . \J : t k ]- Set e = ^ so that (1 + e)r = j> + a — 1. Using that sup fc>1 (tfc — 
^fe-i) < +°° an d the elementary inequality \u + v\ 1+£ < 2 e (u 1+e + v 1+e ) for u, u > 0, 



fc>i 



l+e 



/■+00 

<cy2^m r(1+£) (^T)] + c / E[^( i+£ )(y] 9 (n)rf 
k>i J ° 
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where 5k = k ( 1+£ ) and g is the nonincreasing function defined by g(u) = k ( 1+£ ) on 
[tk-ljtk)- Thus, we deduce from f|3.25j) that, 



Jtk-1 



tfc 



1+e 



< +00. 



It follows from the Chow Theorem (see e.g. |10| ) that (N n ) a.s. converges toward a finite 
random variable which in turn implies by the Kronecker Lemma that 



1 

n ^— < 



n 



tfe-i 



V r (H kT )-V r (Qdu 



tk-i 



->■ a.s. 



Then, (pHHD will follow from 

1 n rtj, 

sup - V / [V(£ fcT ) - V r (du)] du < +00 a.s. 

In order to prove (|3.29|) . we need to inspect two cases for r: 
Case r > 1. 



(3.29) 



We decompose the increment V r (£fc;r) — V r (£ u ) into elementary increments, 



namely 



N(kT) 

V r {t kT ) - V r (Q = V(£ fcT ) - V(^r) + £ V&J - Vfc^). 

£=jV(w)+l 

Owing to the second order Taylor formula, we have for every i£ {N(u) + 1, . . . , N(kT)}: 



where 0/ E (£ ,£„ ). 



Note that a similar development holds for V r (^kT) ~ V r (£kT)- Now, one checks that 
the fact that V G £Q(R d ) implies that \\D 2 V r \\ < C v V r ~ l and that VV is a Lipschitz 
continuous function with Lipschitz constant [vT]i. Consequently 



where we used in the second inequality the standard control \u + v\ s < 2 s ~ 1 (\u\ s + \v \ s ). 
Then, summing over i and using that (VV, b) < j3 owing to (S a>p )(M'), we deduce that 



V r fe T ) - V r (Q < P(kT -u)+ (VV r (Q, a(QdW v ) 

J u 

,-kT 

+C V / V r ~ (Cu)|C«AfcT - £u| 2 + ICvAfcT - £y 



2r 



7iV(i;)+l 



where u = r^^+i- By (S a , p )(i), we can use Lemma EOT m) with 4> = V a and p = s to 
obtain for every s > 0, 



- < C s V^(4) 7 ^ )+r 
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(3.30) 



Applying successively the above inequality with s = 2 and s = 2r > 2 and using the chain 
rule for conditional expectations show that, 



< CSr, A ||-y||„ (l + Efc-i 
for some real constant Cy^iuii . As a consequence, 

sup — 

n>l « 



kT 



V r+a ~ 1 (Cv)dv 



tk 



V 



r+a— 1 



t k -i 



k=l 



E / E * fc -! t yr (^) - 1 + sup-E E*_i 



tk-l 



>1 ™ 



fc=l 



V 



r+a— 1 



tk-l 



Let £G (0, gjg) (note that 



|-(a-l) 



r+a— 1 



£+2(a-i) > ^ since p > 2 and < a < 1). Hence 



(1 + e)(r + a — 1) < p + a — 1 and by Lemma EOT i) and ()3.25|) . one checks that 



+oo 



E 7a+7 E *fe-i 



fc=i 





[/ 







r+a— 1 



tk-i 



tk I.-K 
fc— 1 



+oo 



< 



k=l 



tk-i 



[^ P+a_1 (y] du < +oo a.s. 



by the first part of the lemma. Then, one derives using a martingale argument based 
on (I3.26p . the Chow Theorem and the Kronecker Lemma that 



1 

sup — 

n>l n 



E E ^ 



k=l 



V 



r+a— 1 



{£u)du 



tk-i 



< +oo. 



Case 



< r < 1. 



In that case, we just use that D 2 V r is bounded so that we just have to 



use (|3.3U|) with s = 2 (since a < p + a — 1). This concludes the proof of (ii). 



(in) The fact that a.s., 



1 



f(Qds 



> v{f) is but the statement of Proposition 12.11 



with continuous time notations. Now, let us show that a.s., for every continuous function 
/ such that / = o(F2+ a - 1 ), 



i n 

^E/(W)^-(/). 



(3.31) 



k=l 



First, taking advantage of (|3.27|) . standard weak convergence arguments based on uniform 
integrability show that it is enough to prove that, a.s., (|3.3ip holds for every bounded 
continuous function /. Then, using that weak convergence on M. d can be characterized 
along a countable subset S of Lipschitz bounded continuous functions /, the problem 
amounts to showing that for every Lipschitz bounded continuous function / : R d — > R 



- E /(£ (k-i)r) "^ +00) "(/) a - s - 



(3.32) 



fc=l 



Owing to (Sij.), our strategy here will be to show that almost any limiting distribution of 
the empirical measures is invariant since it leaves the transition operator P T invariant. As 
a first step, we first derive from a standard martingale argument that 



n 

- E /(C(fc-i)r) " E*-2[/(£(fc-i)T)] a.s. 



(3.33) 



k=2 
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Now, we remark that 



Efc-2[/(£(fc-i)r)] = P T f{S(k-2)T) + RkMt(k-2)T) (3-34) 



with R k (x) = E[f(^ (kT) )-f(X^)], (3.35) 



where £ x >T (fcT) denotes the genuine Euler scheme starting from x with step sequence 7 
defined by f|3.24j) . Since / is bounded Lipschitz, 

Rk(x) < CE[|££' 7 - Xt\]1{\ x \<M} + 2 ll/l|ool{|a;|>M} < C M ^/lN{k) + 2 II / II oo l{|a;|>M} , 



(kT) 



where in the second inequality, we used Lemma [3.2l fii) withp = 1. Thus, since 7jv(fc) fc ^ + °°> 
0, it follows from (13. 34ft that, for every M > 0, 

j n 1 ™ 

limsup- Y] (E fc _ 2 [/(C(fc-i)T)] - [/(C(fc-2)T)]) < Climsup- V)l B ( 0) M)=(f(ft-i)T) a - s - 



X n 

Then, it follows from f|3.33|) and from the a.s. tightness of (— ^0,-1^) n >i ^ na *' a - s -' 

n k=i 

n 1 n 1 

- E (/(£(*-dt) - ^/(e (fc -DT)) = -E (/(^-dt) - ^/(e (fc - 2) T)) +o(-) o. 

fc=l fc=2 

Now, since / and P T / are bounded continuous, it follows that, a.s., for every weak limit 
^oo(w,dx) of the tight sequence {rT 1 Y%=i ^ k -i)T)n>i, ^00 (w,/) = u OQ (uj,P T f) for every 
f £ S. This implies that v^uj^dx) is an invariant distribution for P T and one concludes 
the proof by (S^). □ 



4 Rate of convergence for the martingale component 

This section is devoted to the study of the rate of convergence of the martingale (M n ) 
defined in Lemma 13. 11 The main result of this section is Proposition 14.21 where we obtain 
a CLT for this martingale. On the way to this result, the main difficulty is to study the 
asymptotic behavior of the previsible bracket of this sum of four dependent martingales. 
First, we decompose the martingale increment AM„ as follows: 

AM n = E n [A n+1 + B n ] - E n _i[A„+i + B n ], 

where (A n ) is a (C/ n )-adapted sequence defined for every n > 1 by: 

A n = n _i + 4> n - f F (£u)du = I F T (£~ )du- I f F (£,u)du, 

Jl„-1 J(n-3)T Jln-l 

and B n = AM n ^. Keep in mind that E n _i[-B n ] = 0. In the following lemma, we set 

z k ._ x (kT)£ kT yk > 1; 

where, following the notation introduced in ()2.8|) . X^ T ''^ kT denotes the unique solution 
to dY t = b(Y t )dt + a{Y t )dW^ starting from £ kT . 
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Lemma 4.4. Assume b and a are Lipschitz continuous functions satisfying (S a , P ) with 
an essentially quadratic Lyapunov function V and parameters a G (0, 1] and p > 2. Let 
F : B uc ([0, T], M d ) — > M. denote a functional satisfying (Cp) and (Cp). Then, 



1 n 

- ^E fc _ 1 [(AM fe ) 2 ] - E fe _ 2 [(AM; 



2n n— >+oo „ 

fc j J ► a.s. 



(4.36) 



k=2 



and 



1 £ (E fe _ 2 [(AAf fc ) 2 ] - (E fc _ 2 [(E fe C fc+1 ) 2 ] - E,_ 2 [(E fc _ 1 C fc+1 ) 2 ])) ^±^> a.s. (4.37) 



fc=2 



where C k+ \ = A k+ \ + B k with 



)du 



and 



A k+1 = F T (Z k u + 2 )du- f F {Z k u 
Jo Jt 

B k = ~ [ T (Vg F (Zt 2 )^(Zt 2 )dW^ T ). 
Jo 

Proof. We consider the (£? n _i)-martingale (N n ) defined by: 

n 

N n :=Y,T (Vi[(AM fc ) 2 ] - E fc _ 2 [(AM fc ) 2 ]) . 



k=2 



Let e > 0. Using Jensen's inequality, we have 

^E fc _ 2 |AiV fc | 1+£ < G]Tp^E fc _ 2 |Ai\4| 2 ( 1+£ ) < C ^-^E k _ 2 \A k+1 + B k 



|2(l+e) 



k>2 



k>2 



k>2 



Using successively conditional Burkholder-Davis-Gundy and Jensen inequalities and (Cp) 
we have 



Efc-2 \A k+ i + B k 



12(1+=) < 3 i+2 £ C ( 2 ||F|| oo r) 2 ( 1+e ) + (HFiur) 2 ^) 



+T e / E fc _ 2 
4 



|V5 F (^)| 2(1+e) lk(^) 



|2(l+e) 



du (4.38) 



Now, since Vg F is bounded and ||o"|| 2 < CV a , 



E 



fe-2 



|v ffF (e u )i 2(1+£) ik(e u )ii 2(1+£) 



< CE fc _ 2 [U a ( 1+£ )(^)] 

< CCl + Sfe^d^C^r)) (4.39) 



where G k)P {x) = E[ sup U p (^' 7( )]. By Lemma l3.2f i) applied with (ft = V and p 

te[o,T] 

a(l + e) with ee (0, 2=±). it follows that for every k > 2, 



E k _ 2 \AM k \ 2 ^ < C Fi£ , T U a(1+£) (^-2)T)- 
Then, we deduce from Lemma 13.31 applied with 5^ = k~( 1+£ ^ that 

^E fc _ 2 |AiV fc | 1+£ < T^V a{1+£ \!i(k-2)T) < +00 O.8. 



(4.40) 



fc>2 



fc>2 
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since a(l + e) < p + a—1. Finally, using the Chow theorem, it follows that (N n ) is an a.s. 
convergent martingale and the result follows from the Kronecker lemma. 

(ii) Set C k = A k + J3 fc _i. We have AM k = E k [C k+1 ] - E fc _i[Cfc+i] so that 

E fc _ 2 [(AM fc ) 2 ] = E k _ 2 [(E k C k+1 ) 2 } -E k ^ 2 [(E k ^C k+1 ) 2 }. 
Thus, it is enough to show that 



1 " 

- ^E fc _ 2 [(E fc C fc+1 ) 2 ] - E fc _ 2 [(E fe C fc+1 ) 2 ] 



n— >+oo 



» a.s. 



k=2 



1 n 

and - ^E fc _ 2 [(E fe _iC fc+ i) 2 ] - E fc „ 2 [(E fc _ 1 C fc+1 ) 2 ] ^±^> a.s. 



(4.41) 
(4.42) 



k=2 



Let us focus on (|4.4ip . Set q = Using conditional Holder and Jensen inequalities, we 
obtain: 

\E k _ 2 [(E k C k+1 ) 2 }-E k ^ 2 [(E k C k+1 ) 2 }\ = |E fe _ 2 [E k (C k+l - C k+1 )E k (C k+l + C k+l )] \ 
< E fe _ 2 [(E k \A k+1 - A k+1 \ p )*(E k \C k+1 + C fc+1 | 9 )« 

+ E fc _ 2 \(E k {B k+1 - B k+1 ) 2 )^(E k (C k+1 + C k+1 ) 2 )^ 



< (E fc _ 2 |I fc+1 - ^l fc+ i| p )^ {E k _ 2 \C k+1 + C k+1 [ 
+ (E k _ 2 \B k -B k \ 2 y {E k _ 2 (C k+1 + C k+l ) 2 ) 



(4.43) 
(4.44) 



Let us inspect successively the three terms involved in (|4.44p . 

Now set G p (x) = E[ sup V p (Xf)]. Still using Lemma f3.2f i). we show (like previously 

te[0,T] 

for (|4.40p ) that, for every k > 2 and r > 2, 



E fc - 2 |C fc+1 + C k+1 \ r <C(l + G fc _ 2if (£ (fc _ 2)T ) + G|(£ (fc _ 2)T )J < CV3(£ (fc _ 2)T ). (4.45) 
On the other hand since F and / F are bounded Lipschitz continuous functions, 



E k - 2 [\A k+1 - A k+1 \P] <C\\ AE A ._ 2 
< C | E fc _ 2 



sup ie,-^- 2 i p 

u6[(fc-2)T,(fc+l)T] ~ 



sup |fi>-ft>F 

«S[(fc-2)T,(fe+l)T] 



Then, owing to the Markov property, 



+ E fc . 



sup |& - ^~ 2 | p 

«S[(fc-2)T,(fc+l)T] 



A 1. 



E k _ 2 [\A k+1 - A k+1 \P] < C [(H k - 2! 3T,p(kk-2)T) + £k-2,3T,p{Z(k-2)T)) A l] with, 

H ki T, P {x) = E[ sup \C' jikT) ~ C' 1 ^ | p ] and K kiT , p (x) = E[ sup - 

uS[0,T] ~ ve[0,T] 

where (£j' 7< ' )^>o denotes the Euler scheme of X x with step sequence as defined 

by (|3.24p . Now, using that for every v 6 [0, T], 
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x,7 



(feT) 



±x,~f( k ' 1 " 1 \p 



< 2 p sup IC 7 '^' -C' 7(fcT '| p , 

fG[0,T] 
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it follows from Lemma 13.21^) that 



£ — 1 ap 

U k ,T,p{x) < C-i 2 N[kT) V^{x), 



and by Lemma \i.2\ ii). 



so that, for every M > 0, 



Kk,T, P {x) < C{1 + |x| p ) 7 



JV(fcT)' 



(4.46) 
(4.47) 



{\Z(k-2) T \<M} 



+ C1 



{l£(fc-2)T|>AO- 



(4.48) 



Finally, we have 

E k ^ 2 [\B k - B k \ 2 } =E fe _ 2 



Jlu 



a*Vg F (Zt 2 )\ 2 du 



On the one hand Vg F and a being both Lipschitz continuous and Vg F being bounded, 
we have for every x, y £ 

\a*Vg F (x) - a*Vg F (y)\ 2 < (7(1 + \\a(y)f)\x - y\ 2 . (4.49) 

As a consequence, using Schwarz inequality and Assumption (S a , p )(i), it follows that 

E fc _ 2 [|5 fc - 5 fc | 2 ] 



< C E fc _ 2 



1 + sup V 2a {Z*- 2 ) 

u<={(k-l)T,kT] 



sup ie„ - 

, «e[(fc-i)r,fcT] 



fc-2|4 



Owing to Lemma l3.2f i). it follows that 

Efc-2[|-Bfc - B k \ 2 } < CV a (i( k _2)T){U k -2,2T,A{£-(k-2)T) + ^fc-2,2T,4(^(fc-2)r)) 5 ■ 

and by (gig]) and (jUT]) that, 



Efc_2[|-B fc - B k \ ] < CV"(^_ 2 )t) (^ V / 7Af((fc-2)T)l /a (?(fc-2)T) + lN((k-2)T) (1 + |€(Jb— 2)X I ) 
< C" ^((fc-^U + ^ 2a (£( fc -2)T) + IC( fc -2)T| 2(1+a) ), (4-50) 

where we used in the last inequality that V(x) < C(l + |x| 2 ). On the other hand since 
\a*Vg F (x) - a*Vg F {y)\ 2 < C(\\a(x)\\ 2 + |k(y)|| 2 ), 



we deduce likewise from (S ap )(z) and Lemma I3.2N ) that 

Efe-2|B fe - B k \ 2 < CV a (Z {k _ 2)T ). 



(4.51) 



Thus, plugging the inequalities obtained in (|4.45p . (|4.48|) . (|4.50p and (|4.5ip into (|4.44p 
and (|4.43p yields for every M > 0, 

lA(i-A) 

|Efc-2[(EfcC fc ) ]-E fc _ 2 [(EfcC fc ) ]| < CM7^(( fc 2 _2)T) 1 {|5(fe-2)Tl<^}+^ a (^-2)T) 1 {|? (fe _ 2)T |>M}- 

Since ^N(kT) — )• as A: — )• oo and p > 2, it follows that a.s., for every M > 0, 



-i n 1 n 

limsup-^|E fc _ 2 [(lEfcC' fc ) 2 ] -E fc _ 2 [(E fc C fc ) 2 ]| < Climsup - ^ F a (^_ 2)T )l 



fc=2 



n— >+oo Tl 



{\^-2) T \>M}- 



k=2 
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Since p > 2, there exists e > such that a(l + e) < | + a — 1. Hence, it follows from f|3.2Tj) . 
that 

1 71 

sup-£V a ( 1+e >(£ (fc _ 2)T )<+oo a . s . 

n - fc=2 

Then, we deduce by a standard uniform integrability argument that 

(fc-2)r)l{|£ (fc _ 2)T |>Af} 



1 - 

lim sup lim sup — ) j V a (€(k~?fr) 1 ' < , t/i ~ 



M-»+oo n->+oo W 



fc=2 



This completes the proof of (|4.41|) . The proof of (|4.42|) is similar and the details are left 
to the reader. □ 

Lemma 4.5. Assume that b and a are Lipschitz continuous functions, 
(i) For every k > 1, 

E fe _ 2 [(E fe C fe+1 ) 2 ]-E fe - 2 [(lEfe-iC fc+1 ) 2 ] = *(£ (fc _ 2)T ) where, (4.52) 



(a;) = E 



*2T I 2T 



eULI^ -e az\f. 



with At := /„* (F T (X^.) - / F (XJ)) d« 3 t > 0. 

(m) //(Cp) holds, ^ is a continuous function onM. d . As a consequence, if moreover (Cp), 
(S^,) and (S ajP ) /io/d /or a £ (0, 1] and p > 2, 



1 n f 
- ^E,_ 2 [(AM fc ) 2 ] ^±^> o- 2 = / *(ar)i/(d 

71 fc=2 J 



x) a.s. 



(4.53) 



Proof, (i) Let A be a bounded (or nonnegative) Borel functional defined on C( TDl Tr " 1 



Since pathwise uniqueness holds for SDi? (jl.2p (b and cr being Lipschitz continuous), 
there exists a measurable function h : M d x C(M+,R^) — >■ C(R + ,M d ) such that a.s., for 
every fc > 2, Z fc " 2 = x((*- 2 ) T U(k-2)r = h(£ {k _ 2)T , W« k ~ 2 W) (see e.g. p3], Corollary 
3.23). Then, using that £,(k-2)T is ^fe-2 = -7 r (fe_ 2 )T- measura ble, that the Brownian motion 
]Y((k-2)T) j g i nc | e p enc | en t of F^ k _ 2 )T and that, F k r = ^F(k-2)T V J 7 ^^ > one derives 
that 

E fe (A(X« fc - 2 ) T )<^-^)) = E(A(X« fe - 2 ) T )'€(^2)^)|^ ((fc - 2)T) ) 



E[A(X^ k ~ 2 ^' x ) | 



2T 



|z=£(fc_2)3 



Using again the representation with function /i (or the fact that strong uniqueness implies 

(E(A(X« k -Wn\^ (k ~ 2)T) )) 
V V ' J J x& 

where (Xf ) t>0 xe ^d is the flow of 



weak uniqueness), one observes that the spatial process (E^A{X ({k - 2)T) ' x ) \ T^T* ^ 



has the same distribution as \E\A{X X )\T^ T 
SDE (jl.2p at tme 0. Consequently, 



E 



k-2 



E fc (A(X« fc - 2 ) T ^-^) 



(e(A(X x )\^ 



x =Z(k-2)T 



Similar arguments show that 



E x (E(A(X X ) \Ft 
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-W\ 2 



Thus, it follows from the definition of Ak+i and Bk+i that 

E fc _ 2 [(E fc C fc+1 ) 2 ] - E fc _ 2 [(E fc _ 1 C fe+1 ) 2 ] = \E X (e (A^X*) | F^ T f - E (A ,.( A' ' ) | F T ] 
where 

i>2T pzi pr 

L x (X*)x= / F T (X x u+ )du- / f F {X*)du- / <<7*V^(X*),dW„ 

JO JT JO 

Note that the second expression clearly defines a functional on the canonical space. Now, 



X =i(k-2)T 



»2T 



EjE (A^X*) I Ff T f - E (A^X*) I 7* 



an 1 t-w^2 



E 7 



E (A^X*) I F%) - E (Ai (X s ) I Jf )) 



E, 



(E (^(X*) I - E (A^X*) I 



where Ai(X*) = A^X 21 ) - ^ f F (X*)du. The result follows using that Ff = F s , that 
(f^(a*Vg F (X^), dW u )t>o is a (J^-martingale and that EL4 2T — A T \ Ft] = 0. 

(ii) Let x G R d and set 

.) = E(A* T I 7 2T ) - e(a* I 7 T ) - a*Vg F {XZ)dW u . 

J T 

Let (x n ) be a convergent sequence of M. d to x. Owing to the standard identity a 2 — b 2 = 
(a — 6) (a + 6) and Schwarz's inequality, 

|#(x) - *(x n )| < E[|^(X, .) - V(^n, .)| 2 ]^[|^(X, .) + ll>{Xn, .)\ 2 ]^ 

Let 5 = {x} U {x n , n > 1}. Since F and / F are bounded, 



E[|V(x, .) + V(x n , -)| 2 ] < C 1 + supE 



yes 



2T x 2 

a*v 5F (x^)dw; 



T 



p2T 

C(l + sup/ E[|<t*V 5f (X£)| 2 ]^ 
uSS JT 



<C(l + sup sup E[V a {X v u )} < C(l + su P y a (x n )), 

veS ue[T,2T] n>l 

owing to Lemma l3.2f i). Thus, 

|*(x) - *(x n )| < CE[|^(x, .) - r/.(x n , .)| 2 ]5 

and it follows easily that ^ will be continuous if both x 1— > E x \ M\A® T \ F iT \j (i = 1, 2) 

and x 1— )• E z ^ J^ T o"*V5F(X^')dW / n ^ , are continuous. On the one hand F and f F being 
Lipschitz continuous, elementary computations show that for i = 1,2, 

riT AT 2 



E 



/ E [A* I j; T ] du- E [A£ I Ji r ] du 
Jo Jo 



< 



sup \Xf-Xf" 
ie[0,3T] 
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Now, since b and a are Lipschitz continuous functions, for every p > 0, there exists a real 
constant Cfe i0 - iPj r > such that (see e.g. [15] or [29!), 



sup \X?-X?" 
te[o,3T] 



< C b „„ T \x - x n \ 



The continuity of x H > K x (K(A X \ T iT I I , i=l,2, follows. On the other hand using (|4.49|) . 



E 



2T 



2T 



(X7g F (X*),a(X*)dW u ) / (Vg F (X^),a(X^)dW u ) 



T 



,-2T 

\ E\a*Vg F a(X x ) - a*Vg F (X^)\ 2 2 du 
Jt 



<C E[l+ sup |X^| 4 ]5E[ sup \X*-X* n ^]2\<C(l + \x\ 2 )\x-x n 

\ m£[T,2T] «e[T,2T] J 

This concludes the proof. 

Proposition 4.2. Suppose that assumptions of Theorem \2.lV a) hold. Then, 



□ 



^^M(0,a 2 F ) as n^+oo. 



(4.54) 



Proof. By Lemma |4.5 



fc=2 



Then, we only need to prove a Lindeberg type condition (see [10], Corollary 3.1). To be 
precise, we will show that for every e > 



1 n 

-^E^IAM^l^ 



|>e-v/n 



k=l 



n— >+oo 



->• 0. 



First, a martingale argument similar to that of the beginning of the proof of Lemma |3. 11 
yields that 

1 n 

-J2 (E k ^[\AM k \\ lAMkl ^ } ] -E fc _ 2 [|AM fc | 2 l {|AMfc| > £ ^ } ]) 0. 



k=2 



Second, using conditional Holder and Chebyschev inequalities, we have for every e, 5 > 

E,_ 2 [|AM fc | 2 l {|AMfc| > £v ^ } ] < ^^E fc _ 2 [|AM fc | 2 ( 1+<5 )], 
and thanks to (|4.38p and ([4.39p . we deduce that 

E fc _ 2 |AM fc | 2 ( 1+<5 )l {|AA4| > £v ^ } < CG k _ 2Ml+s) (^ k _ 2)T ) < ^V< 1+5 \t {k _ 2)T ). 

£_1 

Thus, taking 5 G (0, 2— ) so that a(l + 5) < p/2 + a — 1, we have for every 5 > 0, a.s., 

n n 
limsup-^E^xOAMfcl 2 !^^,^] < C £ limsup £ V< 1+ V (£ (fc _ a)T ) = 0, 



fc=l 



fc=l 



by applying Lemma EOT m), 



□ 
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5 Study of (B n ,i), (6„ j2 ) and (9 n)3 ) 



In this section, we focus on the remainder terms of the decomposition of the error (see 
Lemma l3.ip . Owing to Proposition 14.21 it is now enough to prove that 



e, 



= — > as n — > +oo for % = 1, 2, 3. 



where — > denotes the convergence in probability. For i = 1,2, these properties are stated 
in Lemma 15.61 and 15.71 For i = 3, the result is obvious. 

Lemma 5.6. Assume b and a are Lipschitz continuous functions such that (S a)P ) with 
parameters a € (0, 1], p > 2, and an essentially quadratic Lyapunov function V satisfying 
liminf| :I .|^ +00 F p+a " 1 (x)/|x| > 0. Let F : H\ C (R + , R d ) -> R be Lipschitz continuous. If 
the step condition \2.11$ holds then 

@n,l P „ 

— > as n — >■ +oo. 



n 



Proof. Since F is Lipschitz continuous, it follows from Lemma \3.2( ii) (applied with p = 1) 
that, for every u£ Ik, 



E 



fc-i 



F T (^)-F T (X^) 



< [i^LipE 



fc-i 



sup iCt ~X t 

te[o,T] 



< Cb,a,T,F ,/7iV(«)(l +^k-l\(,u\ 



Consequently, 



n „ 

\®n,l\ < C bt<TtT ,F^2 I .hN(u)+ldu(l +E fc _i sup |£„|) 
k=l J h V v£l k 
n „ 

< Cb,a,T,F ^2 / HN(u)+ldu (1 + |C(fe-l)Tl) 

fc=1 J/ fe V 

where in the second inequality, we used Lemma [3T2T i). Since liminf| x .|_ s>+00 V p+a ~ 1 (x)/\x\ > 
and N(u) = N(u), we deduce that 



Thus, owing to the Kronecker Lemma 

@n,l n-H 

— — t- u a.s. li 

^ *=! 

with 



a.s. if ^4(l + yP +a - 1 (e (fc _i ) T)) <+oo a.s. 



1 



5fc = Tk \k V^^ du 

Now, as (5k) is nonincreasing, it follows from Lemma l3.3( i) that it is now enough to show 
that Y2k>l ^ k < +°°- We have 



(r+oo 
1+ / 
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7tvh+i 



and 



r [2^n iu < E ^ r-' > du < E ^2m. 

Using that the step sequence (7 n ) is nonincreasing, we deduce from Condition (|2.14|) that 

3 

-du < -%= < +00. □ (5.56) 



Lemma 5.7. Assume b and a are Lipschitz continuous functions satisfying (S aiP ) with an 
essentially quadratic Lyapunov function V and parameters a G (0, 1] and p > 2. Let F : 
B uc ([0, T], M. d ) — > R be a functional satisfying (C F ) and (C F ). If the step condition \2.1J$ 
holds, then 

©n,2 P n 

— — > U as ra — > +00. 
\/n 

Proof. Owing to the Ito formula, we have: 

9 F (£kT) - 9 F (€(k-i)T) = / Ag F (£ u ,£u)du+ {Vg F (£ u ),a(£u)dW u ) where, 

Ag F (x,y) = (Vg F (x),b(y)) + ^Tr (a*(y)D 2 g F (x)a(y)) . 
Then, it follows from the definition of Q n ,2 that 

n r-nT / \ 

©n,2 = ^5 f (6t) - 9 F {Z,{k-i)T) + - Ag F {£ u ,£,u)J du 

r-nT 

+ / (^^u)Vg F ^ u ) - a*^ u )Vg F (^ u ),dW u ). 
Jo 

Since <7 F is bounded, 

g F {inT) -g F (Co) a 



1 

—j= V] g F (t 

v k=i 



fcTj - 5 F (4(fc-i)T) = — 7=-^ ► as n -)• +00. 

vfi 



Then, it is now enough to show that 

f -nT 



Iff - \ L 1 

—= J \Ag F (£u) ~ Ag F (Cu,^u) J du — > as n ->■ +00 (5.57) 



and that, 

4= r T (^*(e«)V^(e„)-ff*(e«)V^(^),£flV u ) ^+0 as n^+00. (5.58) 
V n jO 

First, using that g F is a bounded C 2 -function with bounded Lipschitz continuous deriva- 
tives, that b and a are Lipschitz continuous functions, one checks that 

Ag F {x) — Ag F (x,x) < C (\x — x\.\b(x)\ + \x — x\ + \x — x\ 2 + ||<t(x)|| 2 .|x — x\j . 

Then, using that 

maxM£„-£„|,|£„-£„| J < 2 sup \£v - £,u\, 
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it follows from Lemma l3.2( m) applied with <p = V a , p = 1 and p = 2, that, 

< c {^^vHq(i + 1 6(4) i + ii 2 ) +77v (u)+ i^ a (e«: 

By Asssumption (S a , P ), we deduce that 
1 



E 



nT 



du 


c f 









nT 



e y 2 



a (Cu)] v / lN(u)+ldu. 



Now, since p > 2, |a < p + a — 1, and by (|5.56p . we have 



^ M <cyi< + oo. 



fc>i 



Then (|5.57p follows from Lemma I3.3N ) and the Kronecker Lemma like in the proof of 
Lemma 15.61 

Second, we focus on (|5,58p . Set Zq = and 

n 

Zn = Y J ^= (<r*(yVg F (£ u ) - a*{t u )Vg F {Z u ),dW u ), n > 1. 
fc =i vk J h 

The sequence (^ n ) being a (£/ n )-adapted martingale, it follows from Doob's convergence 
Theorem for L 2 -bounded martingales that (|5.58[) holds if 



supE[(Z n ) 2 ] < +oo. 

n>l 



(5.59) 



Let us show (I5.59p . First, 



k>l 



°*(QVgAtiu)-cT*VgMu: 



du. 



By similar arguments as for (|4.49j) . 

<y*{x)Vg F {x) -a*Vg F (x) 2 < C (l + \\a* (x)f) [\x - x\ 2 + \x - x\ 2 ' 
Then, owing to the fact that u£ \u,u], it follows from (S a , p ) and Lemma l3~2l m) that 

2" 



E 



v*(QVg F (Z u )-a*Vg F (U, 



< C lN{u)+l E[V 2a ^u)}. 



Thus, since u < kT for every u G Ik, 

52E[\AZ H \*]<EUZ 1 \*] + C'E [ E[V 2a (^}^^du 

fc>l fc>2^ 7 fe U 

<c(i + ^ + °°E[y 2a (U]^±i c 



Finally, by a similar argument to (|5.56p . we have 



+oo 



lN(u)+l 



du < C < +°° 



u 



fc>l 



and (|5.59p follows from Lemma [3.3f i) and from the fact that 2a < p+a—1 whenp > 2. □ 
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6 Proof of the main theorems 



The first step for the proof of these theorems is now to state our main result about the 
sequence {V^ n,T \ui, F T )) n >i studied in the two previous sections: 

Proposition 6.3. Let T > 0, a G (0, 1] and p > 2. Assume b and a are Lipschitz 
continuous functions satisfying (S ajP ) with an essentially quadratic Lyapunov function 
V such that liminf N ^ +00 V p+a ~ 1 (x)/\x\ > 0. Assume (S£). Let F : B uc ([0, T], R d ) be 
a functional satisfying (Cp) and (C F ). Finally, assume that the step sequence (7 n )n>i 
satisfies and {2.1$ . Then, 



nT (v^ T \co,F T 



n— >+oo 



(6.60) 



Proof. Owing respectively to Lemma 15.61 Lemma 15.71 and the fact that F is bounded, 
@n.,i> @n.,2 and n +i 3 defined in Lemma 1331 satisfy: 



Qn.,1 + Qn,2 + Qn+1,3 

\J nT 



— > as n — > +oo. 



Then, the proposition follows from Proposition 14.21 and from the decomposition of V^ n ' T ^ — 
¥ U (F T ) stated in Lemma 13. 11 □ 

We are now able to prove Theorems 12.11 and 12.21 

Proof of Theorems 12.1( a) and 12.21 First, let (t k )k>i denote a sequence of positive real 
numbers such that t k — > +oo. Set n k = |_t^J ■ Since F T is a bounded functional, we have: 

TkQ^ J*" F T (^)du - Vu(F T j) - v^T (v^' T \u;,F T ) - F U (F T 



<2||F T || 0O (V^-V^ r ) + ll^r 



tk — nkT k^+c 



■> a.s. 



T DO / rri 

Thus, Theorem 12.21 follows taking u = u. For Theorem I2.1f a), setting ii = uV (\u/T\T), 
and t n = T n , we obtain that 



^(kl " F T(^ V(LfJT)) )^-^(^)) ^AA(0,^)0 as 



n — > +oo. 



Now, 



vl»)(t( u ),F T )-- 



F T (^(lf\T)) )du 



n JO 



N{k) 



< \\F T 



Lr„J 

" lloo \ ^ 

1 n 



fc=l 



and the fact that 7jv(fc)+i < C implies that, 

i=jV(fc-l)+l 

Lr 



7= X] 77v(fc)+i < -7= Yl 

J- r>. , v J- n ■ 



7i- 



fc=i 



i=l 



By (|2.14p and the Kronecker Lemma, 

1 n 

77f = 5Z 7i 



n— >+oo 



> Vs > 3/2. 
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Applying this identity with s = 2 yields the result. □ 
Proof of Theorem 12.1T b). Owing to Theorem 12, Vi a), it is now enough to show that 

Since F T is a Lipschitz bounded functional, it follows from the definition of the previous 
occupation measures that 

\,,e \ pW(X(u;),F T ) - p(- n \Z(cj),F T )} <I^JMp F" E [ sup \^ +s - ^ +s \]du. 

L J Vin JO s£[0,T] 

By Lemma l3.2N t>) and Jensen's inequality, for every q > 1, 

i 

E[ sup |£ M+S - £u +s \\Fu\ < E[ sup |4 +s - Cu+sH-Fu] 9 
se[o,T] se[o,T] 

< c (v*(&4m+i) ' ^ cyf (^ 7 y +1 . 

Thus, we deduce that 

/•r„ 3_i 

/ E[sup |^ +s -4 + s|]^<C^ 7fc 2 *E[W(£ )]. 
JO se[o,T] ^ 

Let 5 be a positive number such that (|2.17p holds. Taking q such that 1/q < S, we deduce 
from (|2.17p and Lemma l3.3( i) that 



3 1 



2 q /"-(-OO 'V ^ 9 

^ VU Jo \/J-iv(u)+i 



dn < +oo. 



We again deduce the result from Kronecker's Lemma. □ 

Proof of Theorem 12.31 We only give the main ideas of the proof of this result about 
the "perfect Euler scheme" (Xt), that is naturally simpler than that of the discretized 
processes. First, the reader can check that setting 

nT 

T— 1 / -iW It 1 •» 1 

u, 



1 r nl 



one obtains a similar decomposition as that of Lemma I3.ll replacing u by u and <p F by <f>p 
defined by 

F (1) = O and / F T {X ( - u) )du ifk>2. 
Jh-i 

The main difference in this decomposition is that the term corresponding to n> i is null. 
Then, since the assumption liminfu.!^^^ 

yp+a i( x )/| x | > o is only needed in the proof 
of the result about 6 n)1 (see Lemma |5T6|) . we deduce that it is not necessary here. Then, 
the sequel of the proof works since the statements of Lemma I3.3I still hold if one replaces 
£ by X. To be precise, the first statements of (i) and (ii) can be directly derived from [25] 
(Chapter 1) and the second ones from an adaptation of the proof of this lemma. □ 
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7 Numerical Test on Barrier Options in the Heston model 



As shown in [22J, our algorithm can be successfully implemented for pricing path-dependent 
options in stochastic volatility models when the volatility process evolves in its stationary 
regime. Furthermore, such stationary versions of stochastic volatility models are more per- 
forming to take into account the behaviour of implicit volatility for short maturities. Then, 
even if the assumptions of our main theorems are usually not satisfied for the functionals 
involved in this context, we choose in this section to illustrate them by such an example. 
To be precise, we test numerically the asymptotic normality obtained in the main results 
on the computation of several Barrier options in a Heston stationary stochastic volatility 
model. The dynamics of the traded asset price process (St)t>o is given by: 



dS t = S t (rdt + V(l - TWtdWl + p^TtdWf), S = s > 0, 
dv t = k(0 - v t )dt + Sy/v~ t dW? , v > 0, 

where r denotes the interest rate, (W , W 2 ) is a standard two-dimensional Brownian mo- 
tion, p € [—1,1] and k, 9 and ? are some nonnegative numbers. This model was introduced 
by Heston ([II]). The equation for (vt) has a unique (strong) pathwise continuous solution 
living in R+. If moreover, 2k9 > q 2 then, (vt) is a positive process (see [32]). In this case, 
the volatility process (vt) has a unique invariant probability uq with gamma distribution, 
namely v§ = 7(0, b) with a = (2k)/^ 2 and b = (2k9)/q 2 . Thus, we assume that (vt) evolves 
in its stationary regime, i.e. that 

£(uq) = fo- 
under this assumption, we showed in [22] that any option premium can be expressed as 
the expectation of a functional of a two-dimensional stationary stochastic process. Let us 
recall the idea: we will write (St) as a functional of a stationary process. Elementary Ito 
calculus yields 

S t = s exp ( r t~lf o v s ds + p J ^T s dW 2 + V±-P 2 j yfcdW}) ■ (7.61) 
Introducing the 2-dimensional SDE, 

\dy t = -y t dt + y/vidWl, 



dv t = k(9 - v t )dt + q^vidWf, 



(7-62) 



and using the fact that 



* 1 f l /"* 9 v t — vq — k9t + k fi v s ds 

y/vldW^ =y t -y + / y s ds and / ^ s dW 2 = ^ , 

Jo Jo S 

we deduce that we can construct a (continuous) map $ from C(M + ,M 2 ) to C(R+,R) such 
that (St)t>o = &((yt,vt)t>o)- Now, we have built (yt) so that (yt,vt)t>o has a stationary 
regime. Denoting by p the invariant distribution of (yt,vt), we obtain that 

E[F(S u 0< t<T)]= E M [Fo <S>((y t ,v t ),0 < t < T)]. 

For further details we refer to |22j . Here, we are interested with an Up-and-Out barrier 
option whose discounted payoff is given by: 

F(S t ,Q< t<T) = e~ rT (S T - K) + l {suPo<t<T s t < L} 
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where L > K > 0. We now specify the discretization. First, the genuine Euler scheme of 
the so-called Heston volatility process (also known as the Cox-Ingersoll-Ross process) (vt) 
cannot be implemented since it does note preserve the positivity. Thus, we must replace it 
by a specific discretization scheme: we denote by (vt) the stepwise constant Euler scheme 
built as follows: 

^r n+1 = K„ + k ln+i {0 ~ v Fn ) + ? v /^(W r 2 n+i - Wl n ) | and v = x > 0. 

Note that convergence properties of this scheme have been studied in a constant step 
framework in [5j (see also [8], PQ, and [2] for other specific discretization schemes). 

Second, we denote by (&) the continuous discretization scheme of (log(||))i>o defined by 
£o = and 

6 = Zv n Hr-±v r jt+ P ^(wi-wl)+J(i - p i)v Tn (wt-w l rn ), t€ [r n ,r n+1 ], n > o. 

(7.63) 

Note that we do not need to introduce the Euler of (yt) since its use is nothing but 
a theoretical way to justify why an algorithm for the approximation of the stationary 
regime can be adapted to this context. Finally, in order to compute the supremum of (£t), 
let us recall the principle of the so-called Brownian Bridge method (transposed to this 
framework). Set 

W t Vn) = P( W Lt ~ W rJ + V^?{Wl +t - Wl) 

and let (Y^ 1 ) denote the Brownian Bridge on [0, 7] defined by Y^' 7 = Wt — ^Wj, t G [0, 7]. 
For every t £ [T n , r n+ i], we have 

6 = e r „ + l r " +i ~S t - r B ) + v^y t w(Tn) ^ +1 - 



Using the independence and the Gaussian properties of the Brownian motion, one deduces 
that, for every n > 1, the processes (£t)te rr^iV-i]' ' G {0, ...,n — 1} are conditionally 
independent given the cr-field c((£ 7; , v r{ , < I < n) and that 

^((6)te[T,^, +1 ]|(^r,>€r I+lJ «r,) = (^Z,^+l,^)) 

= C (x, + ^^*+ ^' 7i+1 , [0, 7;+ j) 

where TV denotes a standard Brownian motion. Then, using the symmetry principle, one 
can show that, for every x,y € M, for every z > max(x,y) and positive A and 7, 

t 2 

P( sup x + (y - x)— + AY t ' 7 < z) = 1 - exp( -^(z - x)(z - y)). 

te[0, 7 ] 7 ' 7^ 

It follows that given (£ r; , £ r(+1 , v rf ), sup t6 [p Jirj+1 ] & can be simulated by the method of 
inversion of the distribution function. 

Let us now detail the algorithm. 

Step 1 : From n = to n = N(T). At each step between n = and n = N(T)—1, simulate 
recursively, v r i and £ r ^ . Then, use the Brownian Bridge method to simulate V n = 
su Pte[r n ,rv +1 ] & given (£ r „ > £r n+1 > ^r„ )• Compute recursively M n := max(Vi, . . . , V n ) = 
max(M„_i, V n ). At time N(T), compute 

= e" rT (s exp(^ T ) - K)+l {so suPtg[0 T] ea . p (6)<L}- 
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Step i : From n = N(T + IVi) + 1 to n = N(T + T;). If M N{T+Ti _ 1+1) = V t -i, re- 
place M N{T+r ._ 1+l) by max(F ij ...,V Ar(T+ri _ 1+1) ). Store , . . . , £ rjV(T+r ._ i)+1 ) and 

(ft, . . . , ^(r+ri-i+i))- As in Ste P !> from n = iV(T + r<_i) + 1 to n = N(T + A), 
compute recursively v T , £ r , and the maximum of Vi, Vi+%, . . . , V n . Then, at time 
N(T + Ti), ' 

+ "T ±i e ~ rT ( s oexp(£ T -Cr.^) -K)+l { sup ao exp(£t-e r )<£} ~ F ) 

1 y te l r i-i' r iV(T+r i _ 1 )+i] 

For the following choices of parameters, 

s = 50, r = 0.05, T=l, p = 0.5, (9 = 0.01, s = 0.1, fc = 2, if = 50, L = 55, 

(7.64) 

we want now to obtain an approximation of the distribution of the (asymptotically normal) 
normalized error 

£n ■= (%((( W ),F)-e- rT E[(5 T -if) + l {suP[) ^ TSl < i} ]) 

First, we need to have an accurate approximation of the (risk-neutral) price. In this way, 
we choose to combine a very long simulation with a variance reduction method taking the 
corresponding Barrier option in the Black-Scholes model as a control variable. Indeed, 
on the one hand, it is well-known that the price of such Barrier option has a closed form 
in the Black-Scholes model (based on the Black-Scholes formula for European options) 
and on the other hand, this price can be approximated using the algorithm described 
above by simply replacing the stochastic volatility (vt) by a constant volatility denoted by 
a. Note that the natural choice for a is the long term volatility which is the mean of 
the stationary volatility process (vt) as well. Then, denoting by (^ 5 ) the genuine Euler 
discretization scheme of the Black-Scholes model (especially with the same trajectory for 
W 1 ) with constant volatility 8, we approximate the price of the option by 

i/M F) - (£ B5 M , F) + (r, Vo, T, K, L) 

where denotes the (explicit) price of the up-and-out barrier option in the Black- 

Scholes model. Doing so with a simulation size N = 2.10 8 , we get the following accurate 
approximation of the premium: 

e- rT E[(S T - K) + 1 {SUP0 ^ T St <L}] « 1, 689. 

Then, setting N = 5.10 5 , we proceed M = 10 4 independent Monte Carlo simulations of Bjy. 
We denote by d 2 the empirical variance of the sample (£^, . . . , £jf ) (which corresponds to 
an estimation of a 2 F ) . In Figure [T] are depicted the density of a centered Gaussian random 
variable with variance a 2 and the empirical density fh (smoothed by a convolution with 
a Gaussian kernel) defined by: 




As a conclusion, this numerical experiment first illustrates that the CLT occurs at a 
reasonable range (for numerical purpose) and also suggests that a local version holds true 
as well ("convergence of the density"). Another extension of our result could be, in the 
spirit of Bhattacharia's result in [3] to establish an invariance principle of Donsker type. 
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Figure 1: Comparaison of the approximate density fh of £^ (dotted line) with the density 
of Af(0,d%), N = 5.10 5 , M = 5.10 3 , h = M~*. 



Appendix 

A. Proof of identity (12351) . 

We have to deduce (|2.19p from f|2. 16[) . First, we have (dropping x in Af): 



E[A 2T \F 2T ] - j 2 J o*Vg F {XZ)dw}j E[A T \F, 



e„ |Ef.i„.i^„i- / rr'Tii. (xv;:u/ir„ ) eu, \f,:\ = e u [e[a 2T \f 2T ]e[a t \f t ]} 

= E U [(E[A T \F T ]) 2] 
since one easily checks that EL4 2T — Aj,|.F T ] = 0. It follows that 

i-2T 



Ta z p = E p 



pZl 

E[A 2T \F 2T \ - J a*Vg F (X u )dW u 



E U [(E[A T \T T }Y}. 



Second, using the Markov property (or the fact that X^' x = cp(X^,W^)) and the 
stationarity of the process, one observes that ELA 2T — A^J-^,] and ELA T |.F T ] have the 
same distribution under ¥ u . In particular, 



(E[A 2T - A T \F 2T }) 



E u 



(E[A T \T T ]Y 



Since E[A, T [J 7 .^] = A T + E[A, T — A^J 7 .^] and E^[^4 T E[^4 2T — A^J^]] = 0, we obtain 
that 



TV; ;=E U [Ai] -2E U 



E[A 2T \T 2T ] £ T a*Vg F (X u )dW u +E U (£ 



, 2T - 2 ' 

a*Vg F (X u )dW u 



All we have to do now is to check that the three above terms correspond respectively to 
the three parts of (|2.19p . First, by Fubini's Theorem, 



E v [A 



T pT 



E U [(F(X^ - f F (X u ))(F(X^ - f F (X v )))]dv. 



>u=0 Jv=0 

Owing to the stationarity of the process under P„, we have 

E V [(F(XM - f F ( Xu ))(F(X^ - f F (X v )))] =C F {\u- v\), 
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JO 



where C F is defined by (I2,20p . This yields 



Second, setting 

we have 

E. 



(7.65) 



C F (u — v)dv du = 2 / (T — u)C F {u)du. 
Jo 

M T = f* <r*Vf(X u )dW u , 
Jo 

r2T -i r2T 

E[A 2T \F 2T ] a*Vg F (X u )dW u = E U [(F T (X^) - f F (X u ))(M?f - M 9 T F )]du 

JT \ JO 

= f T E u [(F T (X^) - f F {X u )){M 9 T l u - M°?)]du 
Jo 

r2T 

+ J E V [(F T (XM) - f F {X u ))(M 9 2 £ - M°r)\du 

Now, the fact that M^+ u — M 9p = g F (XT+ u ) — 9 f {Xt) — Jt +U AgF{X v )dv implies that 
we can make use of the stationarity property to obtain for every u £ [0, T], 

E V [(F T (XW) - f F {X u )){M 9 T l u - M°?)\ 

rT 



= E u (F T (X) - f F (X )) [g F (X T ) - g F (X ) - J Ag F {X v )dv 

= E V [{F T (X) - U{X Q )){M a T * - M g T 1 u )]. 
With similar arguments, one checks that for every u € [T, 2T], 

E U [(F T (X^) - f F (X u ))(M 9 2 f - M^)} = E U [(F T (X) - f F (X ))M g J_ u ]. 
It follows that 



rZl 

E[A 2T \F 2T ] a*Vg F (X u )dW u 

(F T (X) - f F (X )) (tM 9 t p - J* M^_ u du + jf" Ml 



Finally, E v 



g F 

2T-u 



du 



2T ~ 

a*Vg F (X u )dW u 



r2T 




= / E„ 


\a*Vg F (X u )\ 2 


Jt 





TE V [F T (X)Mjf ] 



du 



= T j \a*Vg F (x)\ 2 v(dx) 
owing to the stationarity of the process. This concludes the proof. 



□ 



B. Computation of a 2 when F(a) 



As mentioned in (|2.3p . when <fi = Ah + C, the CLT for marginal functions combined with 
a change of variable yields a 2 = L d |cr*V h(x)\ 2 u(dx) . Let us check this formula starting 
from (|2.16p . Following the notation introduced in (|7.65p . we have 

r2T 

E[A 2T \F 2T ] -E[A T \F T ] - a*Vg F (X*)dW u = ( Pl (x,.)-ip 2 (x,.) where, 



ipi(x, .) 



rT r2T r2T 

/ <f>(X^ +T )du + E[<t>(XZ +T )F 2T ]du- f F (X*)du-(M 9 2 r -M*T) and 
Jo Jt Jt 



E [<P(X* +T ) | Ft] du. 
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In this case, f F = P T (p and using that A and P T commute, one checks that f F — v{f F ) = 
AP T h. This implies that g F = P T <p. For the sake of simplicity we may assume w.l.g. 
v(f F ) = v{<j)) = 0. Then, on the one hand 

ipi{x,.) = / Ah(X u )du + / P u _ T <t>(X 2T )du - g F (X 2T ) - g F (X T ) 

JT JT L 

= h(X 2T ) - h(X T ) - (M 2T - My) + / AP u h(X 2T )du - \g F (X 2T ) - g F (X T ) 

Jo 1 

= h{X 2T ) - h(X T ) - (M 2T - M%) + P T h{X 2T ) -h{X 2T ) - [g F {X 2T ) - g F (X T 

=9 F (X2T) 

= g F (X T ) - h(X T ) - (M 2 h T - Mj). 
On the other hand 

f T P u <t>(X T )du = [ T AP u h(X T )du = P T h(X T ) - h(X T ) = g F (X T ) - h(X T ). 
Jo Jo 



so that 

a 2 F = =M V [(M^ - M^) 2 ] = - / E„[\a*Vh(X u )\ 2 ]du = / \a*Vh(x)\ 2 v(dx). □ 
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